Reentrant behavior of the phase stiffness in Josephson junction arrays 
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The phase diagram of a 2D Josephson junction array with large substrate resistance, described 
by a quantum XY model, is studied by means of Fourier path-integral Monte Carlo. A genuine 
Berezinskii-Kosterlitz-Thouless transition is found up to a threshold value g* of the quantum cou- 
pling, beyond which no phase coherence is established. Slightly below g* the phase stiffness shows 
a reentrant behavior with temperature, in connection with a low-temperature disappearance of the 
superconducting phase, driven by strong nonlinear quantum fluctuations. 
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Two-dimensional Josephson junction arrays (JJA) are 
among the best experimental realizations of a model be- 
longing to the XY universality class and offer the possi- 
bility of controlling and studying a variety of phenomena 
related to both the dynamics and the thermodynamics 
of vortices. In these systems a Berezinskii-Kosterlitz- 
Thouless (BKT) transition separates the supercon- 
ducting (SC) and the normal (N) state, the latter dis- 
playing no phase coherence Q. For a nanoscale size of 
the junctions a new interesting feature shows up in the 
JJA, namely the quantum fluctuations of the supercon- 
ducting phases. These are caused by the non-negligible 
energy cost of charge transfer between SC islands, a con- 
sequence of the small capacitances involved and the fact 
that phase and charge are canonically conjugated vari- 
ables. A relevant effect is the progressive reduction of 
the SC-N transition temperature. Recently, fabricated 
arrays of nanosized junctions, both unshunted Q and 
shunted have given the opportunity to experimen- 
tally approach the quantum (zero temperature) phase 
transition. 

However the mechanism of suppression of the BKT 
in the neighborhood of the quantum critical point and 
its connection with the observed reentrance of the ar- 
ray resistance as function of the temperature is not yet 
clear 0,0,01 In this letter we study the SC-N phase dia- 
gram by means of path-integral Monte Carlo (PIMC) @ 
simulations focusing the attention on the region of strong 
quantum fluctuations, in order to investigate their role in 
suppressing the BKT transition. 

We describe the JJA on the square lattice by a quan- 
tum XY model with the following action 
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where (fij — (pi — tfj is the phase difference between the 
ith and the jth neighboring superconducting islands. We 
assume the presence of weak Ohmic dissipation due to 
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FIG. 1: Phase diagram for a JJA on a square lattice in the 
limit R s > Rq = h/(2e) 2 with r\ = 10~ 2 . The symbols 
are our PIMC results; the solid line reports the semiclassical 
result of the PQSCHA 8]. The dashed line is a guide for the 
eye (see text). 



small currents flowing to the substrate or through shunt 
resistances Q , which reflects into the prescription to con- 
sider the phase as an extended variable 0- The capaci- 
tance matrix reads Cij = C Wj Sij + (z Sij — J^d^i-j+d)] > 
where Co = rj C |jj and C are, respectively, the self- and 
mutual capacitances of the islands, and d runs over the 
vector displacements of the z = 4 nearest-neighbors. 

The quantum dynamics of this system is ruled by the 
competition between the Coulomb interaction of Cooper 
pairs, described by the kinetic term, and the Joseph- 
son coupling represented by the cosine term. The quan- 
tum fluctuations are therefore ruled by the quantum cou- 

(2e) 2 /2C 
It 



pling parameter g = y/ E c /E } , where E c 
is the characteristic charging energy (for rj <C 1) 
is also convenient to use the dimensionless temperature 
t = k B T I E , . In our model Eq. we did not include any 
explicit dissipation term, as dissipative effects are negligi- 
ble provided that the shunt resistance R s 3> R Q g 2 / (27rf), 
where R Q = h/(2e) 2 is the quantum resistance; for 
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smaller R s the Caldeira-Leggett term can be added to 
the action 0, , resulting in a decrease of quantum 
fluctuations. 

Fig. ^ displays our resulting phase diagram together 
with the semiclassical results valid at low coupling |8| . At 
high temperature, the system is in the N state with ex- 
ponentially decaying phase correlations, (t/?i(/?j), and van- 
ishing phase stiffness. By lowering t, for g < g*~3.4, 
the system undergoes a BKT phase transition at t BKT (g) 
to a SC state with power-law decaying phase correlations 
and finite stiffness. When g is small enough (semiclassi- 
cal regime), the critical temperature smoothly decreases 
by increasing g and it is in remarkable agreement with 
the predictions of the pure-quantum self-consistent har- 
monic approximation (PQSCHA) 8]. For larger g (but 
still g < g*) the semiclassical treatment becomes less ac- 
curate and the curve t BKT (g) shows a steeper reduction, 
but the SC-N transition still obeys the standard BKT 
scaling behavior. Finally, for g > g* a strong quantum 
coupling regime with no sign of a SC phase is found. Sur- 
prisingly, the BKT critical temperature does not scale 
down to zero by increasing g (i.e., t BKT (g*) ^ 0): by re- 
ducing the temperature in the region 3.2 < g < g*, phase 
coherence is first established, as a result of the quenching 
of thermal fluctuations, and then destroyed again due to 
a dramatic enhancement of quantum fluctuations near 
t = 0. This is evidenced by a reentrant behavior of the 
stiffness of the system, which vanishes at low and high t 
and it is finite at intermediate temperatures. The open 
symbols in Fig. ^ mark the transition between the finite 
and zero stiffness region when t is lowered. 

These results are obtained using P1MC simulations on 
Lx L lattices (up to L = 96) with periodic boundary con- 
ditions. Thermodynamic averages are obtained by MC 
sampling of the partition function after discretization of 
the Euclidean time u £ [0,/3fiJ in P slices ^f, where P 
is the Trotter number, using the standard Metropolis al- 
gorithm. The actual sampling is made on imaginary- 
time Fourier transformed variables using the same algo- 
rithm developed in Ref. 0: thanks to the fact that the 
move amplitudes are independently chosen and dynam- 
ically adjusted for each Fourier component, this proce- 
dure ensures to efficiently reproduce the strong quantum 
fluctuations of the paths in the region of high quantum 
coupling g. Indeed, test simulations with the standard 
PIMC algorithm showed serious problems of ergodicity, 
though eventually giving the same results. The auto- 
correlation times has been reduced by an over-relaxation 
algorithm [t| over the zero-frequency mode. 

A very sensitive method to determine the critical tem- 
perature is provided by the scaling law of the helicity 
modulus T, a quantity proportional to the phase stiff- 
ness. T measures the response of the system to the ap- 
plication of a twist fco to the boundary conditions along 
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FIG. 2: Size scaling of the helicity modulus Tl at the tran- 
sition temperature. Symbols are PIMC data and the dashed- 
lines are the one-parameter fit with Eq. ©. Left panel: g = 
and t = 0.892 [L Q = 0.456(6)]; right panel: g = 3.4 and 
t = 0.25 [L = 3.32(3)]. The insets show A(t) for different 
temperatures, using the two-parameter fit (see text). 
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where F is the free energy. By derivation of the dis- 
cretized path-integral expression of the partition func- 
tion, the PIMC estimator for T is easily obtained, in 
analogy to that of Ref. ^(J- Kosterlitz's renormalization 
group equations provide the critical scaling law for the 
finite-size helicity modulus T^: 
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where Lq is a non-universal constant. Following Ref. 
the critical temperature can be found by fitting T l (t) [t 
vs L for several temperatures according to Eq. © with 
a further multiplicative fitting parameter A(t). In this 
way, the critical point can be determined by searching 
the temperature such that A(t BKT ) — I, as illustrated in 
FigEl This is the technique we used to get the filled sym- 
bols in Fig. ^ Using this procedure the critical tempera- 
ture can be determined with excellent precision. For in- 
stance in the classical case we get £ BKT (.g=0) = 0.892(2), 
in very good agreement with the most accurate results 
from classical simulations Also in the regime of 

strong quantum coupling, g — 3.4, the PIMC data for 
Tz,(i BKT = 0.25) are very well fitted by Eq. 0, as shown 
in Fig. [3] Moreover, this figure points out the sensi- 
tivity of this method to identify t BKT : at temperature 
higher (lower) than the critical one the helicity modulus 
decreases (increases) much faster with L than T^(i BKT ). 
At higher values of the quantum coupling, g > g* , the 
helicity modulus scales to zero with L — > oo and P — > oo 
at any temperature. 
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FIG. 3: Helicity modulus T_l divided by the best fit with the 
expression © for g = 3A and different temperatures: A, O; 
and □ correspond to t — 0.2, 0.25, 0.3, respectively. 
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FIG. 5: Trotter-number extrapolation of Tl for g — 3.4. 
Two series of data for t = 0.1 (•) and 0.2 (■) are reported, 
for four different lattice sizes: from the top to the bottom 
L = 8, 10, 12, 14. The lines are weighted quadratic fits. 




FIG. 4: Temperature behavior of the helicity modulus Tg(t) 
on a 8 x 8 lattice, for different values of g. The data are results 
from the Trotter extrapolation. 



FIG. 6: Finite-size scaling of the helicity modulus Tl for 
g = 3.4 at fixed P — 101. The lines are guides for the eye. 



At variance with the standard BKT theory, in the 
regime of strong quantum fluctuations we find a range 
of coupling values, 3.2 < g < g*, in which the helicity 
modulus displays a non-monotonic temperature behav- 
ior. In Fig. 01 Ti(i) is plotted for different values of g 
on the 8x8 cluster: up to g = 3.0 it shows a monotonic 
behavior similar to the classical case, where thermal fluc- 
tuations drive the suppression of the phase stiffness. In 
contrast, for g~3A, the helicity modulus is suppressed 
at low temperature, then it increases up to t ~ 0.2; for 
further increasing temperature it recovers the classical- 
like behavior and a standard BKT transition can still 
be located at t ~ 0.25 (Figs. El and EJ. A reentrance 
of the phase stiffness was found for a related model in 
Ref. [lfj, but the authors concluded that the low tem- 
perature drop of the helicity modulus was probably due 
to having a finite Trotter number. In order to ascertain 
this point, we have performed systematic extrapolations 
in the Trotter number and in the lattice size, as illus- 
trated in Figs. El and El for g = 3.4. We have found 
no sign of anomalies in the finite-P behavior: the ex- 
trapolations in the Trotter number appear to be well- 



behaved and already in the expected asymptotic regime 
0{1/P 2 ) 0, for P > 60 (Fig. El). The extrapolation to 
infinite lattice-size shown in Fig. El clearly indicates that 
Tl scales to zero at t — 0.1, while it remains finite and 
sizeable at t = 0.2. Hence, the outcome of our analysis 
is opposite to that of Ref. [10(, i.e., we conclude that the 
reentrant behavior of the helicity modulus is a genuine 
effect present in our model, rather than a finite- Trotter 
or finite-size artifact. 

In order to understand the physical reasons of the reen- 
trance observed in the phase stiffness, we have studied 
the following two quantities: 

( cos (pij(u)) , (4) 
K = {(fij( u )-^ij) 2 ) . (5) 

with ifij — Jq 13 du (fiij(u) and ij nearest-neighbor 

sites. The first quantity is a measure of the total (thermal 
and quantum) short-range fluctuations of the Joseph- 
son phase: in particular, it has a maximum were the 
overall fluctuation effect is weakest. The second quan- 
tity represents instead the "pure quantum" spread of the 
phase difference between two neighboring islands and has 
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FIG. 7: Top panels: vs t; bottom panels: (cos tpij(u)} 
vs t. The quantum coupling is g — 1.0 in the left panels and 
g = 3.4 in the right panels. The circles are PIMC data and 
the dashed lines are PQSCHA results. 



been recently studied in the single junction problem [14j; 
more precisely, measures the fluctuations around the 
"static" value (i.e., the zero- frequency component of the 
Euclidean path) , it is maximum at t = and tends to zero 
in the classical limit (i.e., g/t — > 0). 

The quantities and JSJ on a 8 x 8 lattice are com- 
pared in Fig. 0for two values of the quantum coupling, 
in the semiclassical (g — 1.0) and in the extreme quan- 
tum (g = 3.4) regime. In the first case (cosipij(u)) 
decreases monotonically by increasing t and the pure- 
quantum phase spread A^, shows a semiclassical linear 
behavior which is correctly described by the PQSCHA. 
At variance with this, at g = 3.4, where the reentrance 
of T(i) is observed, (cos (pij(u)) displays a pronounced 
maximum at finite temperature. Besides the qualitative 
agreement with the mean- field prediction of Ref. (lif . 
we find a much stronger enhancement of the maximum 
above the t = value. This remarkable finite-i effect can 
be explained by looking at A^(g=3.4) (Fig. top-right 
panel) . Its value is an order of magnitude higher than the 
one in the semiclassical approximation and, notably, it is 
strongly suppressed by temperature in a qualitatively dif- 
ferent way from A^(g-l.O): the pure-quantum contribu- 
tion to the phase fluctuations measured by A^, decreases 
much faster than the linearly rising classical (thermal) 
one. Thus the interplay between strong quantum cou- 
pling and temperature turns out in a finite-i minimum 
of the total fluctuations of the Josephson phase. This 
single-junction effect in a definite interval of the quan- 
tum coupling (3.2 < g < 3.4) is so effective to drive the 
reentrance of the phase stiffness. 

As for this low-temperature transition, the open sym- 
bols in Fig. represent the approximate location of the 
points (t,g) where T(t) becomes zero within the error 
bars: in their neighborhood we did not find any BKT-likc 



scaling law. This fact opens two possible interpretations: 
(i) the transition does not belong to the XY universality 
class; (ii) it does, and in this case the control parameter 
is not the (renormalized) temperature, but a more invo- 
lute function of both t and g. Further investigations are 
needed to answer this question. 

In summary, we have studied a model for a JJA in the 
quantum fluctuation dominated regime. The BKT phase 
transition has been followed increasing the quantum cou- 
pling j up to a critical value g* ~ 3.4 where t BKT ~ 0.25; 
above g* no traces of BKT critical behavior have been 
observed. Remarkably, in the regime of strong quantum 
coupling (3.2 < g < 3.4) phase coherence is established 
only in a finite range of temperatures, disappearing at 
higher T, with a genuine BKT transition to the normal 
state, and at lower T, due to a nonlinear quantum me- 
chanical mechanism. 
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